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PACS 74 . 72 . -h - Cuprate superconductors 

PACS 74 . 20 . Mn - Nonconventional mechanisms 

PACS 79.60.-i - Photoemission and photoelectron spectra 

Abstract. - Recent angle resolved photoemission [1] and scanning tunneling microscopy [2] mea- 
surements on underdoped cuprates have yielded new spectroscopic information on quasiparticles in 
the pseudogap phase. New features of the normal state such as particle-hole asymmetry, maxima 
in the energy dispersion and accompanying drops in the spectral weight of quasiparticles agree 
with the ansatz of Yang et al. for the single particle propagator in the pseudogap phase. The 
coherent quasiparticle dispersion and reduced asymmetry in the tunneling density of states in the 
superconducting state can also be described by this propagator. 



The anomalous electronic properties of the pseudogap 
phase in underdoped cuprates continue to challenge the- 
ory. Angle resolved photoemission (ARPES) measure- 
ments of the Fermi surface in the pseudogap state show 
only disconnected Fermi arcs centered on the nodal di- 
rections. [3-5] A recent high resolution ARPES study on 
underdoped BSCCO by Yang and collaborators [1] exam- 
ined the spectra below and above the chemical potential fi 
by carefully dividing out the Fermi function. They found 
that the spectra displayed particle-hole asymmetry away 
from the nodal direction in the pseudogap state but not 
in optimal doped samples, so that the asymmetry is a 
property of the pseudogap phase. Scanning tunneling mi- 
croscopy (STM) by Kohsaka and collaborators [2] has de- 
termined the coherent quasiparticle dispersion in a wide 
range of hole densities in the superconducting state, by 
an ingenious analysis of the voltage dependent spatial in- 
terference patterns. [6] The asymmetry in the tunneling 
density of states measured by STM is weaker due to the 
inherent particle-hole symmetry of superconductivity. Re- 
cent angle integrated photoemission (AIPES) experiments 
on underdoped samples [7] found the density of states var- 
ied linearly with energy in contrast to the constant value 
that appears at a band edge in two dimensions. In this let- 
ter we analyze ARPES, STM and AIPES results using an 
ansatz for the single particle propagator proposed earlier 



by three of us. Key properties of the quasiparticles (QP) in 
this propagator, such as maxima in the energy dispersion 
in the normal pseudogap state, the connection between 
the QP dispersion in the superconducting state and the 
hole density, are directly confirmed by these experiments. 
An alternative explanation of particle-hole asymmetry in 
ARPES experiments has been proposed by Anderson [8]. 

As the hole density, x, crosses the underdoped pseu- 
dogap region of the phase diagram, the cuprates crossover 
from a full Fermi surface metal at overdoping to a Mott in- 
sulator at stoichiometry. A Mott insulating state is driven 
by strong onsite Coulomb repulsion and not per se by 
translational symmetry breaking. The stoichiometric un- 
derdoped cuprate YBa 2 Cu40s, shows no sign of a charge 
or spin density modulation on cither Cu or O sites. [9] It 
is therefore desirable to examine theories where the trun- 
cated Fermi surface appears as a precursor to the Mott 
state without translational symmetry breaking. An ex- 
ample is the 2D array of 2-leg Hubbard ladders studied 
by Konik, Rice and Tsvelik (KRT). [10] They showed that 
lines of zeros in the propagator G R (k, uj = 0, x) associated 
with the charge gap at half-filling, enclosing a commensu- 
rate area of 1 el/site, do not move with light doping. The 
Luttinger Sum Rule on the area enclosed by the condition 
Re[G R (k,uj — 0,a;)] > is fulfilled with the Fermi surface 
truncated to small hole pockets. 
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Table 1: The bare hopping parameters t' , tg and the superex- 
change J are quoted in units of the nn value to = 0.3eV. The 
hole density determines the pocket area and is consistent with 
the value determined independently [1]. 
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A relevant question is whether similar behavior can oc- 
cur in a fully 2D system. Honerkamp and coworkers [11] 
argued for such behavior from the similarities in the flow of 
the response functions in functional renormalization group 
analyses of Hubbard models near half-filling. Yang, Rice 
and Zhang [12] (YRZ) proposed an adaptation of the KRT 
approach to the lightly doped t-J model. They introduced 
a self-energy, H R (k,u>,x), which diverges at u = on a 
surface spanned by elastic particle-particle umklapp scat- 
tering analogous to the behavior of the ladder model. In 
the 2D square lattice this umklapp surface is a diamond 
connecting antinodal points on the Brillouin zone bound- 
ary. Note this umklapp surface appears as the energy gap 
surface also in the case of wider Hubbard ladders with 
more than 2 legs. [13] 

YRZ took over the form for the gap function A R (k, x) 
from the renormalized mean field theory of Zhang et al. 
[14] for the resonant valence bond (RVB) state [15, 16] of 
the strong coupling 2D t-J model and proposed an ansatz 
for the single particle propagator: 



G R (k, lu, x) = g\x)/[u - & - w, x)] + G, 



(1) 



where Y, R (k,u, x) = \A R (k,x) 
self-energy. The energy £^ 



\ 2 /(uj + ^ 0) ) is the RVB 
= — 2t(x)(cos k x + cos k y ) 



vanishes on the umklapp surface. The renormalized dis- 
persion £k = — At'(x) cos k x cos k y — 2i"(x)(cos 2k x + 
cos2k y ) — (j. includes hoppings out to 3 rd nearest neighbor 
with hole density dependent coefficients t(x) = g t (x)to + 
3g s (x)Jx/8, t'(x) = g t {x)t' , and t"(x) = g t (x)t^. The 
Gutzwiller factors g t (x) = 2x/ (1 + x) and g s (x) = 
4/(1 + x) 2 account for the no-double occupation projec- 
tion on the kinetic and superexchange terms in the t-J 
model. The RVB energy gap takes the form A R (k,x) = 
A(°)(x)(cos k x — cos ky). \ corresponds to the weakly dop- 

>. Pa- 



ing dependent homogeneous amplitude < c\„^,cr 
rameters are quoted in Table 1. 

As discussed by YRZ the Fermi surface in the pseudogap 
regime with divergent G R (k,oj = 0,x) has 4 pockets cen- 
tered on the nodal directions enclosing a total area related 
to the hole density. [12] ARPES experiments [1,3-5] show 
arcs rather than closed pockets, but as YRZ pointed out 
the QP weight in G R varies rapidly with very small values 
on the back sides of the pockets, which can account for the 
failure to observe these parts of the pockets. Many experi- 
mental features, e.g. the slow variation of the nodal Fermi 
velocity vp and Drude weight scaling with hole density, are 
reproduced by the YRZ propagator. [12] The umklapp sur- 
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Fig. 1: (Color online) Comparisons between (a) QP dispersion 
Ek, (b) spectral weight Zh, from the YRZ propagator and the 
values obtained from the ARPES results by Yang et al. [1] 
Error bars reflect the uncertainties in the fitting procedures; 
Zk at low energies in cuts 2, 3 have large error bars due to 
uncertainties in the choice of the rising background. 



face along which the energy gap opens up, lies above /i, 
causing particle-hole asymmetry in the QP spectra in the 
pseudogap state. The spectra are symmetric only along 
the nodal directions where A R (k,x) — 0, and the asym- 
metry increases away from the nodal directions. 

The spectra observed by Yang et al [1] along three cuts 
(see inset in Fig.l) in the pseudogap phase (T=140K) 
are shown in the Brillouin zone in Fig. 2(a). They are 
fits to a broadened Lorentzian plus a background term 
to be defined below. The fitting uses an implementation 
of the nonlinear least-square Marquardt-Levenberg algo- 
rithm. We determine the QP dispersion, Ek, weight Zk, 
and inverse lifetime 7& for each k- value from the form. 



A(k,u) 



Im[z k /(Lu - E k + i~f k )}+ A B (k,u>) (2) 



For the background, momentum distribution curves 
(MDC) of cut 1 are well fitted with an almost con- 
stant background A B (k,ui), for energy distribution curves 
(EDC) of cuts 2 and 3 we use a broadened step function 
which rises at lower energies well below /z, 



A B (k,uj) 



y k (tanh [A fc (O fc - lo)] + 1) + a k (3) 



The fitted MDC(EDCs) for the three cuts are shown in 
the upper panels of Fig. 2(a), and the two components, 
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Fig. 2: (Color online) Fits to the ARPES spectra along the 
cuts (1-3) (see inset of Fig.l) [1]. In (la) the MDC of cuts 1 
and in (2/3a) the EDC of cuts 2 and 3 are shown, the exper- 
imental data are the solid black lines, and the fitted A(k,u) 
are the dashed red lines. Individual components of the fits are 
displayed in (b) for QP Lorentzian peaks and (c) for smooth 
background A B (k,u>). 



the Lorentzian peak and background, in the lower panels 
(b/c). For cut 1 along the nodal directions we note that 
the fitting to a broadened Lorentzian with an almost con- 
stant background works well for all energies. For cuts 2-3 
EDC fitting, as fc moves away from hp, the background 
component A B (k,oj) increases at lower energies. We be- 
lieve this increase arises from contamination from adjacent 
fc-values, due to the enhanced broadening at low energies. 
Since our interest focusses on energies near /i we believe 
this choice of background which is constant near mu is 
reasonable. 

Fig.l shows the results for the QP properties of inter- 
est. The QP dispersion, Ek in Fig. 1(a) rises linearly with 
k along the nodal cut 1. Particle-hole asymmetry is evi- 
dent in cuts 2 (3) with a maxima in Ek lying above (be- 
low) /i, respectively. Note that these maxima are not at 
the boundary of the reduced Brillouin zone as the case 
would be in the presence of a broken translation symme- 
try, a point emphasized by Yang et al. [1] The QP weight, 
Zfc, represented in the dispersion plot as the line width, is 
plotted in Fig. 1(b). We see an almost constant Zk for cut 
1 along the nodal direction, but a fall off at the maximum 
in Ek for cuts 2 and 3 in agreement with the YRZ ansatz. 
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Fig. 3: (Color online) The red curves are the contours of the 
constant QP energy below and close to the Fermi energy. The 
black dots are the turning points deduced from STM interfer- 
ence data. The dashed blue lines are the hole Fermi pocket in 
the normal pseudogap state. The cyan curves are the turning 
points of the iso-energy contour of the corresponding band in 
the SC state. Note only the turning point of the iso-energy 
contour with the largest spectral weight are shown. 



The inverse lifetime 7& is essentially constant indicating 
substantial inelastic scattering at this elevated tempera- 
ture of 140K. To summarize with reasonable parameters 
the phenomenological YRZ propagator gives a good fit to 
the asymmetry of ARPES spectra and describes both the 
maxima in the QP dispersion, Ek, away from the nodal 
direction and the accompanying drop off in the weight, Zk- 
In their original discussion of the ARPES spectra, Yang 
et al. [1] presented results for a fourth scan closer to the 
antinodal point (cut 2 in Fig. 4 of Ref. [1]) than scans (1-3) 
in Fig. 3 of Ref. [1] . This showed particle-hole symmetry 
in the accessible energy window near the Fermi energy 
(< 40 meV). Note, this energy window is too small to see 
the actual gap on the particle side. However, the obser- 
vation that the minimum in the gap is located right at 
the chemical potential led Yang et al. [1] to suggest that 
the pseudogap at the antinodal regions has a pairing ori- 
gin, in agreement with the proposals of Kanigel et al. [4] 
and Anderson [8]. In a standard pairing theory the en- 
ergy gap opens up along the Fermi surface whereas in the 
YRZ propagator the gap opens up at a fixed surface in k- 
space. As such some of use (YRZ) argue that this occurs 
due to the presence of umklapp particle-particle scatter- 
ing in addition to the Cooper channel, similar to the case 
of the half-filled Hubbard 2-leg ladder. [11] As a conse- 
quence of this difference, the RVB energy gap, A R , acts as 
a charge gap and only the quasiparticle excitations at the 
Fermi pockets contribute to the reduced Drude weight and 



p-3 



Kai-Yu Yangl,3, H.-B. Yang2, P. D. Johnson2, T. M. Ricel,3, Fu-Chun Zhang3 



c/3 
O 

n 




0.30 
0.25 
0.20 



g 0.15 



T3 

01 
(0 

O 

% 0.10 
C 

0.05 

cs 

g 0.00 
-0.05 



-0.10 



"5 arc-tip^ (theory) 1 

■ hole pocket area (theory) 

V Bogoliubov arc (expt) 

y the simple Luttinger theorem (expt) 



v 



8 
▼ 



0.00 0.05 0.10 0.15 0.20 

doping x 



0.25 



Fig. 4: (Color online) The density of states in the super- 
conducting state for various dopings. The d-wave SC gap is 
chosen to saturate at ends of Fermi pockets at a fixed value 
As| max = 0.14 to agree with experiments. Two choices for the 
RVB gap A are displayed (red A = 0.2, blue A = 0.5). 
The dashed black curves show the STM data [2] in units of to 
(300meV). The inverse of lifetime T = a\E QP \ with a = 0.1 is 
used. 



reduced coherent stiffness in the superconducting state. 
That are characteristic of the pseudogap state. [17] 

At this point we wish to comment the apparent contra- 
diction between the hole Fermi pockets associated with the 
YRZ form and the electron Fermi pockets deduced from 
recent quantum oscillation experiments on underdoped 
YBCO samples. [18]. These experiments are performed 
at high magnetic fields and low temperatures whereas 
the ARPES experiments of Yang et al. [1] were made at 
zero field and at a relatively high temperature of 140K. 
Measurements of the Hall constant on the YBCO sam- 
ples show a hole-like behavior under the latter conditions 
but electron-like Hall constant under the conditions of the 
quantum oscillation experiments [19]. This suggests that 
the sign change in the Hall constant and the electron pock- 
ets, arise from a reconstruction of the Fermi surface [20] 
at high magnetic fields, possibly related to the similar sign 
change in the Hall constant that accompanies the super- 
lattice charge and spin order in the static stripe phase. [19] 

The recent STM results of Kohsaka et al., [2] requires 
a generalization of the YRZ propagator to the supercon- 
ducting state. Originally YRZ added a self-energy term 
in G R (k, u, x). [12] However, during this work we realized 
that this is appropriate only in the presence of particle- 
hole symmetry. It is more convenient to rewrite Eq.l 



G R (k,Lu,x) = J2 WS/(w - ££) + Gi 



Fig. 5: (Color online): Hole densities implied by areas en- 
closed by different interpretations of the QP dispersion shown 
in Fig. 3, and dispersion extrapolated using LDA band struc- 
ture (Yfrom Ref [2]), hole pockets determined by turning points 
in Fig. 3 and lines connecting antinodal points (\y from Ref [2], 
□ from this work), cyan line and ■ from the YRZ propagator 
for the Fermi pocket. . 
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where E k = £fc ± 

&4 0) + |A2f and (W^)- 1 = 1 
The generalization to the superconducting state is 
straightforward by treating the system as a two- 
band d-wave superconductor with gap Ag(fe, x) = 
As{x)[cos(k x ) — cos(k y )] (for simplicity we use a single 
SC gap) leading to a propagator 

G s (k,Lo,x)=J2wZ/[Lo-E%-\A s (k,x)\ 2 /( L o + E%)] (5) 



(4) 



with four QPs E iyk = ± V / (E%) 2 + A s (fc, x) 2 with spectral 
weight z it k- 

The superconducting gap opens up along the hole pock- 
ets. Kohsaka et al [2] derived the coherent QP dispersion 
in the superconducting state from the evolution of the in- 
terference pattern in the STM spectra measured in a wide 
surface area as the tunneling voltage is changed in BSCCO 
samples, ranging from near optimal doping (x — 0.19) 
to strongly underdoped [x — 0.06). This determines a 
set of (/-vectors connecting the turning points in the iso- 
energy contours of the QP dispersion. [6,21] These turn- 
ing points can be compared to the values obtained from 
Eq.5. We use the same values of the parameters to, t' Q , 
t'o and J as fits to the ARPES spectra at x = 0.12. Fol- 
lowing Zhang et al [14] we assume the RVB gap value 
A°(x)[— A°(l — x/x c )] drops linearly with x with a crit- 
ical doping x c = 0.20. Note the value A = 0.5 is used 
in fitting to ARPES results on a sample with x = 0.12. 
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In fitting to the STM results we shall allow for the possi- 
bility that the superconducting As and RVB A gaps are 
not independent and that A may change in the super- 
conducting state from its value in the normal state. 

The remaining issue is the form for the superconduct- 
ing gaps Ag(k, x). Generally we expect that Ag(fe, x) will 
not continue to grow as k moves beyond the Fermi pock- 
ets and the normal state QPs drop below the Fermi en- 
ergy. Therefore we assume a d-wave form for As(fc, x) = 
A°(x)(cos k x — cos k y ) only for fe-values on the pockets and 
set As(k,x) to be constant beyond the pockets. This re- 
sults in a two-gap description [22-30] with the antinodal 
energy gap dominated by the RVB gap at underdoping 
while the superconducting gap along the hole pockets is 
more important near optimal doping. We display in Fig. 3 
the evolution of the contours of constant QP energy near 
/i. Note the strong variation in QP weight discussed ear- 
lier, means that only the outer contours closest to the 
zone center are relevant. The turning points of these con- 
tours are determined by the maxima in the inverse velocity 
(dE(k)/dk±)~ 1 along these contours. The comparison to 
the STM results is illustrated in Fig. 3 for 4 representative 
hole densities. We see that the overall agreement is quite 
satisfactory. Note that the contours defined by the turning- 
points are slightly larger than the underlying hole Fermi 
pockets and enclose a slightly larger area. This causes 
the hole density implied by the turning point contour to 
be slightly larger than the actual hole density. This dis- 
crepancy is evident in Fig. 3(a) of Kohsaka et al. [2]. The 
relevant areas are plotted in Fig. 5. 

Note that in the above analysis we used a standard 
d-wave form for the superconducting gap, A$(k,x) = 
A s (x)(cosfc x cos/uy) as observed in ARPES. Kohsaka 
et al. in their STM experiments found a different angular 
dependence with a strongly reduced gap near to the nodal 
directions. The origin of this discrepancy is not clear at 
present. One possibility could be enhanced forward impu- 
rity scattering in the samples used in the STM experiments 
which Haas et al. showed can lead to a gap suppression 
near to the nodal directions. [31] However since both tech- 
niques are sensitive only to the surface layers it is far from 
clear that such an explanation is viable. 

Kohsaka et al. [2] report an abrupt end to the coherent 
QP dispersion in the STM experiments when k reaches 
the umklapp surface. The YRZ propagator contains co- 
herent QP also in the gapped antinodal regions. QP in 
these regions, however, should be sensitive to the strong 
local variations in the hole density and the accompany- 
ing large variations in the RVB energy gap, A^(a;). [32] 
This may well lead to localization of the corresponding 
QP states which lie at energies away from the chemical 
potential fi. The constant value of \i makes the Bogoli- 
ubov QP states associated with the superconducting gap 
much less sensitive to disorder. This difference offers a 
plausible explanation for the change from propagating to 
localized states that Kohsaka et al. [2] propose takes place 
as k reaches the umklapp surface. 




Fig. 6: (Color online): The upper panels are the angle inte- 
grated spectra I(E,T,x) near the chemical potential (at zero 
temperature) for various doping (0.05 - 0.14). The lower pan- 
els are the first order derivative of I(E,T,x), with the peak 
position E p shifting towards negative electron energy (i.e. pos- 
itive hole energy), opposite to the shift in the chemical poten- 
tial. Blue, black, cyan and red curves are for the tempera- 
tures fcsT = 3, 6, 10, 15 meV. The same parameters as shown 
in Table 1 are used, except A = 0.4. Experimental data for 
La2Cu04 can be found in Fig.xxx of Hashimoto et al [7]. 



Kohsaka et al. [2] reported characteristic tunneling 
spectra, which we can compare to the DOS obtained 
from the superconducting propagator in Eq.5 N(u) = 
J2i k z i,k&Ei fc-u»j as shown in Fig. 4. The experimental re- 
sults are shown in Fig. 4 as the dashed black curves. At 
higher hole densities, 0.1 < x < 0.19, the spectra are 
dominated by the maximum superconducting gap which is 
larger than the RVB gap, A°(x). As a result the QP bands 
near the Fermi energy are all split by As(fe) with symmet- 
ric low energy DOS. At lower hole density the RVB gap 
A (a;) rises and exceeds the maximum of Ag(fc) leading to 
two structures in the DOS. One at lower energy related to 
As(k,x) at the tips of arc, and second at higher energies 
associated with the RVB gap. If we keep the RVB gap A 
to the value used in the ARPES fits, the DOS shown in 
blue in Fig. 4 displays a much stronger particle-hole asym- 
metry than the experiments. However, if we reduce the 
value of A to 0.2, the agreement is much improved, sug- 
gesting that the two gaps are not independent of each 
other. Note that the total gap near the antinodal points 
does not change so much. We do not claim a quantitative 
fit to the STM DOS but the main features are reproduced 
at least qualitatively, by the YRZ ansatz. 

Next we turn to recent angle integrated photoemis- 
sion (AIPES) experiments on hole doped La2Cu04 which 
showed several anomalous features at underdoping in the 
normal pseudogap phases. First the total DOS (i.e. an- 
gle integrated) show an approximate linear dependence on 
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Fig. 7: (Color online): The shifting value of the peak position 
E p of the derivative of dl jdE towards positive hole energy for 
varous doping and temperature. It shows a clear quadratic 
temperature dependence. 



the hole energy, E, with only a small quadratic correction 
which grew with increasing x. This behavior contrasts 
with the constant DOS expected from a small hole pocket 
in a valence band in two dimensions. Secondly, the shift 
of the peak in the first derivative of the angle integrated 
spectra (I(E,T,x)), dl/dE oc T 2 /T co h approximately and 
moves to positive hole energies (not negative). The value 
oiT co h(x) increases with increasing, x. Interestingly, these 
anomalous features are reproduced by the YRZ model dis- 
persion for quasiparticles. In Fig. 6 (a) we show the re- 
sults for the angle integrated spectra I(E, T, x) calculated 
by multiplying the total DOS for quasiparticles by the 
Fermi function. An approximate linear dependence on E 
is evident. In Fig. 6 (b) we show that the peak E p of 
the derivative | dl/dE | moves to positive hole energies 
opposite to the negative shift of the chemical potential, 
with a roughly quadratic dependence on the temperature, 
Fig. 7. The resulting characteristic coherence tempera- 
ture T co h{x) increases with increasing x. The key anoma- 
lous properties of the AIPES spectra are well reproduced 
by the YRZ form. 

The YRZ ansatz for the single particle propagator ex- 
plains the anomalous properties of the underdoped pseu- 
dogap phase as a precursor to the Mott insulating state at 
stoichiometry. The recent spectroscopic measurements on 
underdoped BSCCO using ARPES and STM have shed 
new light on the evolution of the QP properties. In this 
paper we have shown that the YRZ ansatz provides the 
basis to understand the key features of the new experimen- 
tal data such as particle-hole asymmetry, energy disper- 
sion and wave vector dependent spectral weight of normal 
state QP as well as the coherent QP dispersion and DOS 
in the superconducting state over a range of hole densities. 
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